Density-induced BCS to Bose-Einstein crossover 
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We investigate the zero-temperature BCS to Bose-Einstein crossover at the mean-field level, 
by driving it with the attractive potential and the particle density. We emphasize specifically the 
role played by the particle density in this crossover. Three different interparticle potentials are 
considered for the continuum model in three spatial dimensions, while both s- and d-wave solutions 
are analyzed for the attractive (extended) Hubbard model on a two-dimensional square lattice. For 
this model the peculiar behavior of the crossover for the d-wave solution is discussed. In particular, 
in the strong-coupling limit when approaching half filling we evidence the occurrence of strong 
correlations among antiparallel-spin fermions belonging to different composite bosons, which give 
rise to a quasi- long- range antiferromagnetic order in this limit. 



I. INTRODUCTION 

The BCS to Bose-Einstein (BE) crossover has been 
widely studied in the last several years, |l]-|(| being orig- 
inally motivated by the occurrence of a short coherence 
length in high-temperature superconductors. The evolu- 
tion from large overlapping Cooper pairs (BCS limit) to 
small nonoverlapping bosons (BE limit) has essentially 
been envisaged by relying on the associated two-body 
problem in the three-dimensional case, wherein bound 
fcrmion pairs (composite bosons) form as soon as the 
strength of the attractive interparticle potential exceeds 
a threshold. The emphasis on the role of the interparti- 
cle potential has, however, somewhat overshadowed the 
effects of the particle density on the crossover itself, even 
though on physical grounds one would expect both the 
interparticle potential and the density to play an essen- 
tial role. The role played by the density is suggested 
especially when one analyzes the experimental phase di- 
agram of the high-temperature cuprate superconductors 
in terms of the BCS-BE crossover, since in this case it 
would be the (effective) carrier density (that is related to 
the doping level) to drive the system from the vicinity of 
the BE (underdoped) to the BCS (overdoped) limit. 0|] 

Purpose of this paper is to study the combined effects 
of the particle density and the interparticle potential on 
the (zero-temperature) BCS-BE crossover, in order to 
characterize how physical quantities evolve by varying, 
in particular, the particle density. To this end, we will 
set up a "phase diagram" in the space of the potential 
strength and of the density, where the locations of the 
alternative BCS-like, crossover-like, and BE-like regions 
will be identified for several types of potentials. 

Previous work on the BCS-BE crossover has utilized: 
(i) a three-dimensional contact potential for the contin- 
uum model (ii) the separable potential introduced 
by Nozieres and Schmitt-Rink (NSR) for the three- 
dimensional continuum model, with a characteristic mo- 
mentum cutoff k Q [^,^); and (Hi) & negative-/7 Hubbard 
model on a two-dimensional square lattice, with cither 



an on-site Q or a nearest-neighbor |9||lC|] attraction. In 
this paper, we shall consider both the continuum model 
in three spatial dimensions using three different types 
of interparticle potentials, and the s- and d-wave solu- 
tions for the attractive (extended) Hubbard model on a 
two-dimensional square lattice. This will enable us to 
study the effects of the particle density on the BCS-BE 
crossover in a rather systematic way. 

For the continuum case, it turns out that the finite 
range of the potential allows for the occurrence of the 
density-induced BCS-BE crossover, which is instead not 
possible in the case of a contact (zero-range) potential. 
As a consequence, the size of the BCS-like region in the 
"phase diagram" gets progressively enlarged by increas- 
ing the range of the interparticle potential. 

For the lattice case, the shape of the "phase diagram" 
and the physical interpretation of the alternative regions 
therein depend markedly on the symmetry (s or d) of the 
pairing. In particular, for the d-wave pairing an increas- 
ingly larger range of correlations among the composite 
bosons sets up when approaching half filling, thus estab- 
lishing a tendency toward the formation of a quasi-long- 
range-ordered antiferromagnetic state. In addition, the 
eZ-wave pairing, being associated with an interaction of 
a finite range on the lattice, enables the density-induced 
BCS-BE crossover to occur over a wider range of the pa- 
rameters with respect to the s-wave pairing, in analogy 
to what found for the continuum case. 

All results presented in this paper have been obtained 
within a zero-temperature broken-symmetry mean-field 
approach, which thus appears capable of producing re- 
sults that are sensible on physical grounds under widely 
different physical conditions. 

The plan of the paper is as follows. We discuss the 
three-dimensional continuum case in Section II and the 
two-dimensional lattice case in Section III. Section IV 
gives our conclusions. 
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II. THREE-DIMENSIONAL CONTINUUM CASE 

In this Section, we examine the three-dimensional 
BCS-BE crossover in the continuum case for three types 
of interaction potentials, and determine how the range 
of the potential influences the coupling vs density "phase 
diagram". Specifically, we consider the contact and the 
separable NSR potentials mentioned in the Introduction, 
plus a non-separable Gaussian potential. 

For the three-dimensional continuum case, a contact 
potential has often been adopted as the reference model 
that captures the essence of the expected physics of the 
BCS-BE crossover, as a function of the coupling strength 
for given particle density. For a contact potential, the an- 
alytic solution at the (zero temperature) mean-field level 
and with the inclusion of Gaussian fluctuations has been 
determined, [Tl| for all values of the coupling strength 
(regularized in terms of the scattering length a of the 
associated two-body problem) and of the density (repre- 
sented in terms of the Fermi wave vector fcp). All relevant 
physical quantities can thus be expressed in terms of the 
dimensionless parameter fcjro, with Uf being positive by 
definition and a changing its sign as soon as a bound 
state develops. For this reason, by keeping kp fixed and 
varying a from — oo to +00 one can pass with continuity 
from the BCS to the BE regime across the crossover re- 
gion; on the contrary, by keeping a fixed and varying kp 
one is not able to pass from the BCS to the BE regime, 
since the parameter kpa cannot change its sign in this 
way. 
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FIG. 1. Phase diagram (g,kpa ) for a contact potential 
in three dimensions, as obtained from the available analytic 
solution (see text). 

Figure 1 shows the (g,kpa Q ) "phase diagram" for 
a contact potential in three dimensions, as obtained 
from the available analytic solution, where BCS-like, 
crossover-like, and BE-like regions are identified by draw- 
ing the two curves corresponding to kp^ pa ir — 27r and 




kF^pair = 1/tt. Here, £ pQ i r represents the correlation 
length for pairs of opposite-spin fermions, while the two 
values (27r, 1/tt) of the parameter kp£ pair characterize, 
in the order, the lower limit of the BCS-like region (with 
large overlapping Cooper pairs) and the upper limit of 
the BE-like region (with small nonoverlapping bosons), 
the crossover region being constrained in between. |l2] ] 
The (effective) coupling constant in this case is repre- 
sented by g — exp (a G /a), where a is an arbitrary unit 
of length. As anticipated above, it is evident from this 
"phase diagram" that for a contact potential in three di- 
mensions it is not possible to cross over from the BCS 
to the BE region by varying the density alone at fixed 
coupling strength. 

In the three-dimensional continuum case, therefore, to 
examine the density-induced BCS-BE crossover a poten- 
tial with finite range in real space is required, or equiv- 
alently, it is necessary to introduce a momentum cutoff 
k Q in momentum space. In this context, one may utilize 
the separable potential V(k, k') = Vw(k)w(k') (between 
fermions with opposite spins) adopted in Ref. Q, with 
V < and w(k) = (1 + ((kl/fco) 2 )^ 1 / 2 , and later used by 
some authors. [|||(| Since this potential may yield unphys- 
ical results [Q (also because the factorization w(k)w(k') 
is somewhat arbitrary), we have considered in addition a 
non-separable potential V(k, k') = V(k — k'), which we 
have taken for definiteness of the Gaussian form: 



V(k-k') = ^exp{-|k-k'| 2 /fco} 



(1) 



(V < 0). By doing so, it will also be possible to de- 
termine how alternative ways of introducing an effective 
range in the interaction potential affect the boundaries 
of the crossover region in the (IVI/V^, kp/k D ) "phase di- 
agram" (V c being the critical value of \V\ for which a 
bound state appears in the two-body problem). In ad- 
dition, and contrary to the separable NSR potential, the 
Gaussian potential (|l|) leads to momentum-decoupling 
effects for small values of k D , which have recently been 
proposed as characteristic features of high-temperature 
superconductors. jl3) 

For the NSR and Gaussian potentials, an analytic so- 
lution for the (zero-temperature) BCS-BE crossover at 
the mean-field level is lacking. For these potentials, we 
have thus solved numerically the coupled equations for 
the gap function A(k) and the chemical potential /x: 



A(k) 



and 



(2tt)3 V[ ' > 2E(k>) 
dk! (. £(k') 



(2tt) 3 



1 - 



E(k')J ' 



(2) 



(3) 



where £(k) = k 2 / (2m) —jj, (m being the fermionic mass), 
E(k) = ^/£(k) 2 + A(k) 2 , and n is the particle density. 
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We recall that, while for a separable potential the gap 
function acquires the form A(k) = A u>(k), for a non- 
separable potential A(k) does not follow in general the 
wave vector dependence of the potential. In addition, for 
the NSR separable potential analytic expressions for A Q 
and \x can be obtained in the two (BCS and BE) limits. 
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FIG. 2. Phase diagram (\V\/V C , k F /k a ) for the (a) NSR 
and (b) Gaussian potentials in three dimensions (see the text 
for the meaning of the different curves). 

Figure 2 shows the (\V\/V c ,kp/k ) "phase diagram" 
for the (a) NSR and (b) Gaussian potentials, where the 
two characteristic curves kp^ pa ir = (27T, 1/vr) for each po- 
tential have been identified as for Fig. 1 . In addition, we 
have reported in Fig. 2 the curve corresponding to /.t = 
(broken line) for both potentials. By comparing Fig. 2 
with Fig. 1, we note that the boundary between the BE- 
like and the crossover-like regions is not much altered by 



the introduction of a finite cutoff fc„; on the contrary, 
the boundary between the BCS-like and the crossover- 
like regions is drastically modified, with the BCS-like re- 
gion extending even to values |V| >> V c for sufficiently 
high densities. Note also the reentrant shape of the curve 
kptipair — 27r at kp ~ k for both potentials, which makes 
it possible to cross over from the BCS to the BE region 
by varying the density for fixed \V\ > V c . 

It is further interesting to note that the three curves 
of Fig. 2a and of Fig. 2b depart from the common origin 
(|y|/I4 = l,kp = 0). Near this origin, in fact, kp « k Q 
and the gap A„ = A(k = 0) on the BCS side is propor- 
tional to (kp/2m) ex-p{ir/ (2k pa)}, as given by the solu- 
tion for a contact potential. PJTl|] Since ^ pa i r oc kp/A a 
in the BCS limit, (l^] keeping the product kp^ pa ir oc 
exp{— -rr/(2kpa)} — constant, requires a — > — oo when 
kp — ► 0. This, in turn, implies that the curves corre- 
sponding to kF^p a i r ~ 27r depart from the point |V| = V c 
on the kp = axis. On the BE side, on the other hand, 
when kp << k Q (and |V| > V c ), £, pa ir coincides with 
the bound-state radius r D ; O] keeping thus constant the 
product kF^pair = kpr when kp — > implies r — > oo, 
i.e., \V\ - V c . @ 

The reentrant shape of the curve kpS,pair — 27r in Fig. 2 
at kp ~ k Q , too, can be understood by simple analytic 
arguments as follows. When kp << k a , the expression 
A Q //i oc exp{7r/(2A:i?a)}, which is valid on the BCS side 
for a contact potential, can again be used. At fixed value 
of |V| < V c (such that the scattering length a is neg- 
ative), A Q //i vanishes when kp — * 0. In this way, one 
approaches the (weak-coupling) BCS limit for decreas- 
ing kp at fixed a < 0. When kp >> k a , on the other 
hand, it is necessary to distinguish the NSR from the 
Gaussian potential. For the NSR potential, the value of 
A(kp) = A a w(kp) can be obtained analytically in the 
BCS limit (and kp >> k Q ), in the form 



A(*j 



exp 



1 



N V(k F ,k F ) 



(4) 



where N a — mkp/(2ir 2 ) is the density of states at 
the Fermi level (per spin component) and V(kp,kp) ~ 
(k /kp) 2 V. In this case, the decrease of \V(kp, kp)\ for 
increasing kp jk Q overcomes in Eq. (^) the increase of N a , 
and drives the system toward the BCS (weak-coupling) 
limit. For the non-separable Gaussian potential, on the 
other hand, it is the effective reduction of the density 
of states and not the decrease of the potential strength 
to drive the system toward the BCS limit for increasing 
kpjk at given V. To verify this statement, we recall 
that in the BCS limit E(k!) in Eq. (0) is strongly peaked 
about kp. When kp >> k a , A(k) is thus also strongly 
peaked about kp owing to the form (|l]) of the potential. 
In this way, for given value of k, the integral over k' in 
Eq. (j^) extends effectively over a sphere centered about k 
with radius of the order k a . The relevant density of states 
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gets thus reduced from the value N Q by a geometrical fac- 
tor R of the order (Awk^/S)/ (8nkpk ), which represents 
the ratio of the effective volume of integration and the 
BCS spherical shell of width 2k . We then obtain in the 
BCS limit for the Gaussian potential when kp » k a : 

A (fc F )cx^exp{-i-j. (5) 
to l_ RJy V J 

In this case, it is thus the decrease of R for increasing 
fei? to drive the system toward the BCS (weak-coupling) 
limit for a given value of V. 

In summary, we have shown that, although the qualita- 
tive behavior of the curves corresponding to kp£,pair = 2tt 
in Figs. 2a and 2b looks similar, the physical mechanism 
behind them appears to be quite different. For the sepa- 
rable NSR potential the increase of kp/k results into a 
reduction of the interaction strength, while for the non- 
separable Gaussian potential it results into a reduction 
of the relevant density of states. These two alternative 
effects allow the density-induced BCS-BE crossover to 
occur in the two cases, respectively. 

The curves of Fig. 2 have been drawn cautiously, up to 
values of \V\/V C for which the condition \fj,\ « kg /(2m) 
is satisfied. This condition avoids, in fact, instabilities of 
the system in the bosonic limit, which unavoidably occur 
when a fermionic potential with finite momentum range 
k D is considered, and are due to the boson-boson effective 
interaction potential having a dominant attractive part 
in this case. [Q,[l5| We have verified, in particular, for 
the Gaussian potential that the bosonic chemical poten- 
tial hb = 2/i + e (where e is the bound-state energy of 
the associated two-body problem) ||[|] becomes negative 
when |//| ~ kg/m, a behavior which can be attributed to 
an overall attraction among the composite bosons (with 
the compressibility being, however, still positive). In this 
context, we mention also that the existence of a compe- 
tition between pair and quartet condensation in a Fermi 
liquid with a finite-range attraction has recently been in- 
vestigated. [[l6| 

A related instability toward phase separation (with 
the compressibility becoming instead negative) has been 
generically pointed out for the attractive (extended) 
Hubbard model on a two-dimensional square lattice, |^] 
which we pass now to examine in detail in the context of 
the BCS-BE crossover. 

III. TWO-DIMENSIONAL LATTICE CASE 

In this Section, we examine the BCS-BE crossover for 
a two-dimensional attractive Hubbard model, again ad- 
dressing specifically the role played by the particle den- 
sity in driving this crossover. To consider the d-wave 
besides the s-wave solution, we take the fermionic po- 
tential to contain an inter-site besides an on-site term. 



In addition, for the d-wave solution we adopt two dif- 
ferent single-particle dispersions, in order to mimic the 
low-energy electronic band structure observed for the 
cuprates in different doping ranges, Jl7],[ll| as discussed in 
more detail below. We recall that the issue of the <i-wave 
symmetry in the context of the BCS-BE crossover has 
been discussed briefly in Ref. Jl9| and more extensively 
in Refs. § and 

For this model we thus take: 

V(k,k') = U + 2V(cos(k x - k' x ) + cos(k y - k' y )) (6) 

where V < 0. In particular, when V = we con- 
sider an on-site attraction U < 0, with the hopping in 
the kinetic term of the fermionic Hamiltonian limited 
to nearest-neighbor sites (which corresponds to the or- 
dinary negative-?/ Hubbard model |5)); when V < we 
consider instead an on-site repulsion U > 0, with the 
hopping in the kinetic term of the fermionic Hamilto- 
nian either limited to nearest-neighbor sites or ranging 
over second- and third-neighbor sites (which corresponds 
to an extended attractive Hubbard model, with an on- 
site repulsion pi]). Recall that the term proportional 
to V in Eq. (||) is associated with an attraction between 
opposite-spin fermions on neighboring sites in the square 
lattice. 

A. s-wave solution 

The negative-?/ Hubbard model (obtained by taking 
V = and U < in Eq. plays on the lattice an 
analogous role to the contact potential in the contin- 
uum case. For this model, the only nontrivial solution 
to the gap equation (|^) (with the integration over the 
wave vector being now limited to the two-dimensional 
Brillouin zone (BZ)) has s-wave symmetry (A(k) = A D ), 
and the dispersion relation £(k) = — 2i(cos k x + cos k y ) + 
nU/2 - fi (t > 0) contains the Hartree shift nU/2. |§ 
The corresponding (\U\/t,n) "phase diagram" is shown 
in Fig. 3a, with the boundaries between the alternative 
(BCS, crossover, and BE) regions identified like for the 
continuum case. Note that the critical value U c for which 
a bound state appears in the two-body (n = 0) problem 
consistently vanishes for the s-wave solution in two di- 
mensions. Note also that the reentrant shape of the curve 
kp£,pair = 27r (as well as of the curve kp^ pa i r = 1/tt) 
stems now from the fact that the van Hove singularity of 
the density of states is approached when n tends to 1 (half 
filling). (Owing to the symmetry of the density of states 
about half filling, the diagram of Fig. 3a is also symmetric 
about half filling.) Note finally that the crossover from 
the BCS to the BE region by varying the density for fixed 
\U\ is possible only for \U\ ~ 2. it, the BE and crossover 
regions being depressed in this case to extremely small 
values of n. 
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FIG. 3. (a) Phase diagram (\U\/t,n) for the s-wave so- 
lution of the negative-?/ Hubbard model in two dimensions 
( with on-site attraction and nearest-neighbor hopping); (b) 
Phase diagram (€ ,k%/(2m)) for a contact potential in two 
dimensions (as obtained from the available analytic solution). 

It is interesting to compare the "phase diagram" of 
Fig. 3a for the two-dimensional negative-?/ Hubbard 
model, with the "phase diagram" (e ,kp/(2m)) (in units 
of (m«^) _1 , where a is again an arbitrary unit of length) 
for the contact potential in the two-dimensional contin- 
uum case, for which the analytic solution is also available. 
|p3| , p| This "phase diagram" is shown in Fig. 3b, where 
now the boundary curves with kp^ pa i r = (2-7T, l/7r) cor- 
respond to straight lines. [We have verified numerically 
that the solutions to the coupled equations (||) and (||) 
for the negative-?/ Hubbard model and for the contact 
potential in two dimensions coincide within a few per- 
cent when n ~ 0.1, even for values of \U\/t of the order 



of some units.] Note that the absence of a threshold for 
the occurrence of a bound state in two dimensions makes 
the density-induced BCS-BE crossover possible even for 
a short-range potential, in contrast to the behavior for a 
three-dimensional contact potential obtained previously 
(cf. Fig. 1). 

A comment on the nature of the bosonic limit for 
large values of \U\/t when approaching half filling in 
Fig. 3a is in order at this point. For the continuum 
case, the limit of a "dilute" gas of composite bosons can 
be reached for any particle density, insofar as the range 
of the residual bosonic interaction (or, equivalently, the 
size of the composite bosons) vanishes for large values 
of the fermionic interaction strength (barring the insta- 
bility problem mentioned at the end of Section II). In 
the lattice case, instead, the lattice spacing provides an 
additional length scale in the problem, which makes it 
possible to depart from the "dilute" gas limit irrespec- 
tive of the size of the composite bosons. It is, in fact, the 
"overlap" of the centers of mass of the composite bosons, 
which is forced by increasing the density on the lattice, 
to make the usual "dilute" gas condition kpa << 1 (as 
defined for the continuum model) no longer represent- 
ing a "dilute" gas situation in the lattice case. For such 
a high-density gas of composite bosons, therefore, the 
underlying fermionic degrees of freedom are expected to 
become significant again. 

To make these arguments more quantitative, let us 
consider the following boson- like operator: 

bt = $>(k)4 T ct_ ki (7) 

k 

where cj_ creates a fermion with wave vector k and spin 
er, <?(k) represents the pair wave function, and the sum 
over k is limited to the Brillouin zone in the lattice case. 
The ensuing commutator 

[bo A] = E l-9( k )| 2 (! - ™ki - n-ki) (8) 

k 

with riko- = c^o-Cko-, can be regarded a c-number provided 
(n\s.a) << 1 over the relevant set of states. [Q The nor- 
malization condition 

E ^ = N ( 9 ) 

k.o- 

(N being the total number of fermions), implies that 
(^k<r) > n/2 over a finite region of k space, whenever 
the range of the sum over k remains limited (as for the 
lattice case). This restriction, in turn, implies that the 
commutator (|^) cannot be considered as a c-number, as 
soon as n is an appreciable fraction of the unity. Note 
that the size of the composite bosons does not enter the 
above argument. 

In summary, we have argued on quite general grounds 
that, to reach a satisfactory bosonic limit in the lattice 
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case, the condition kp^ pa i r < 1/n valid in the continuum 
case has to be supplemented by the condition n << 1. 



B. d-wave solution 



When V < and U > in Eq. (|6|), a d-wave solution 
of the type A (k) = Ai(cosfe x — cosfc y ) can be considered. 
In this case, we take either a nearest-neighbor dispersion 
relation £(k) = -2t(cosfc x + cosk y ) + n(U + W)/2 - fx 
as before (t — V model) , or a second- and third- neighbor 
dispersion relation £(k) = At' cos k x cos k y + 2i"(cos 2k x + 
cos2k y ) + n{U + AV)/2 - fj, with t' > and t" > 
(t' — t" — V model). In both cases we have included 
the Hartree shift. |2^] The latter dispersion relation fa- 
vors the formation of bound pairs with d-wave symme- 
try at low density, pi]] the critical value V c (for a bound 
state to appear in the two-body problem) vanishing when 
t" < 0.5t'. For the former dispersion relation limited to 
nearest-neighbor sites, on the other hand, a finite value 
of V c occurs also for the d-wave solution at zero density. 

i 

As already mentioned, the two different single-particle 
dispersion relations that we have adopted are meant 
to mimic the low-energy electronic band structure for 
the cuprates in different ranges of the doping level 5. 
Specifically, the (t' , £")-dispcrsion captures the small-arc 
features of the Fermi surface detected in underdoped 
cuprates, 0,^5| for which we can interpret n ~ S ~ 0.15. 
The i-dispcrsion, on the other hand, reproduces the main 
features of the large Fermi surface and its doping depen- 
dence for nearly optimally doped cuprates, ] fl8| provided 
one interpretes n ~ 1 — 5 with 0.15 ~ 5 ~ 0.30. 

The "phase diagrams" for the t — V and t' — t" — V 
models are shown, respectively, in Figs. 4a and 4b, where 
the boundaries between different regions have been iden- 
tified as in Fig. 3a. |26| Note in Fig. 4a the occurrence 
of a finite critical value V c (V c /At = 1.83), past which 
the BE region develops for n = (V c vanishes instead 
in Fig. 4b). We also mention that the BCS-like region, 
corresponding to small values of |V| in Figs. 4a and 4b, 
consistently supports nonvanishing values of the gap Ai, 
contrary to a recent statement for the t — V model. |2~ 
Note finally that the density-induced BCS-BE crossover 
is allowed for the t' — t" — V model but not for the t — V 
model. This finding is, in turn, consistent with our previ- 
ous assertions that the t' — t" — V model might be relevant 
to the underdoped range of the cuprate superconductors 
and that the BCS-BE crossover scenario might apply to 
that range. 

A comment on the additional region at the upper-right 
corner of the "phase diagrams" in Figs. 4a and 4b, which 
is delimited by a long-dashed line and has been identified 
as "correlated bosons" (CB), is in order. 
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FIG. 4. Phase diagram (|V|/4t, n) for the d-wave solution 
of the extended Hubbard model with attraction between near- 
est-neighbor sites in two dimensions, considering (a) near- 
est-neighbor or (b) second- and third-neighbor hopping. 



This region is associated with a peculiar behavior of 
ipair for the d-wave solution as a function of |V| when n 
j approaches half filling, in the sense that tl pa ir does not 
show a monotonic decrease for increasing \V\ and con- 
verges asymptotically (when |V| — * oo and n < 1) to a 
finite value, which is larger than the lattice spacing. 

This asymptotic value increases with n and eventu- 
ally diverges as n — » 1, making thus the product kp^ pa i r 
arbitrarily large. [This peculiar behavior is absent for 
the s-wave solution discussed previously, for which ^ pa i r 
is instead a monotonically decreasing function of \V\ 
for any given n. Consistently, the CB region is miss- 
ing in Fig. 3a.] On physical grounds, the divergence of 
ipair cannot be attributed either to the system converg- 



es 



ing to a BCS-like regime with large overlapping Cooper 
pairs or to the size of the composite bosons becom- 
ing infinitely large. Rather, the divergence of ^ pa ir is 
due to the establishing of quasi-long-range-order corre- 
lations among the composite bosons, which reside indi- 
vidually on nearest-neighbor sites. Under these circum- 
stances, £,pair weights preponderantly the correlation be- 
tween opposite-spin fermions belonging to different com- 
posite bosons, rather than the usual intra-boson correla- 
tion. Accordingly, for the d-wave solution when n — ► 1, 
£pair no longer represents the radius of the composite 
bosons in the strong-coupling limit. 

To make the above argument more quantitative, we 
recall the original definition of Cvair in terms of the two- 
particle correlation function: 12 3] 



Jdr gu (r) r 2 
I dr 5T,i( r ) 



(10) 



where 



SUW = ^ <V>J(r)Vl(0)lMOM(r)> 
(#|V}(r)^(0)|$> 2 



n- 
1 



(11) 



is the pair-correlation function for opposite spin fermions 
(the last expression holding specifically for the BCS 
ground state |$)). From this definition, it appears ev- 
ident that, in principle, j.(r) does not distinguish be- 
tween opposite-spin fermions belonging to the same pair 
or to different pairs. In practice, in the strong-coupling 
limit g^i (r) represents (the square of) the pair wave func- 
tion whenever no correlation is established among the 
composite bosons; in this case, £ po ir tends to the bound- 
state radius, as one verifies for the s-wave solution. In 
the case that a definite correlation is established among 
the composite bosons, on the other hand, <?Td( r ) spreads 
over a large (and even infinite) number of lattice sites, 
and S, P air increases (and eventually diverges) accordingly. 
In this case, <?|,;(r) embodies the correlation among dif- 
ferent composite bosons and is totally unrelated to the 
pair wave function. 

The occurrence of this novel feature for the d-wave so- 
lution is evidenced in Fig. 5a, where the amplitude 



Rn) = T7 ex P{-« k ' R n} 



BZ 



A(k) 
2E(k) 



(12) 



is reported for \ V\/t >> 1 and n = 1 over a grid of lattice 
sites R„ (TV being the total number of sites); in this way, 
0(R„)| 2 represents the lattice version of <?t,l( r ) given by 
Eq. (pi]). We see from this figure that the amplitude 
( |l2| ) develops clear structures near the lattice diagonals, 
thus establishing definite correlations among opposite- 
spin fermions belonging to different composite bosons. 
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(b) 



FIG. 5. (a) d-wave amplitude 0(R n ) over the sites R n of a 
two-dimensional square lattice, obtained for the t — V model 
when \V\/t » 1 and n = 1; (b) Pictorial representation of 
the ordering of the composite bosons with d-wave symmetry 
on the two-dimensional square lattice when \V\/t » 1 and 
n = 1. 

This behavior can also be checked analytically, making 
use of the fact that at half filling /i equals the Hartree 
shift (U + 4V)/2 and thus the single-particle disper- 
sion £(k) is negligible in comparison to |A(k)| when 
\V\/t >> 1 (barring values of k for which A(k) vanishes 
by symmetry). In this case 



A(k) 



A(k) 



2E(k) 2|A(k) 



= isgn(cosk x 



cosk y ) (13) 



in Eq. (|12|) alternates sign according to the d-wave sym- 
metry, and its lattice Fourier transform 



(Rn) 



111 

tt2 



-1) 



n x -\-n v 



(14) 
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(with R„ = (n x ,n y ) in units of the lattice spacing) de- 
cays as a power law for increasing distance when n x +n y is 
an odd integer, while it vanishes identically when n x + n y 
is an even integer. [In contrast, for the s-wave solution 
A(k)/£J(k) tends to a constant value when \V\/t >> 1 
and n = 1, and the corresponding Fourier transform is 
a Kronccker delta <Sr„,o-] The ensuing picture bears 
strong resemblance with an antiferromagnetic ordering 
on a square lattice, with opposite spins alternating over 
the two interpenetrating sublattices in which the square 
lattice can be partitioned. The ordering of the composite 
bosons (in the strong-coupling limit when approaching 
half filling), as envisaged from the above considerations, 
is shown schematically in Fig. 5b. Note that the quasi- 
long-range order associated with the algebraic decay in 
Eq. ( |l4|) corresponds to a divergent t; pa ir, even though 
it does not represent a true long-range antiferromagnetic 
order of the spins. Physically, such a strong correlation 
among fermions with opposite spins stems from the orig- 
inal fermionic attraction between opposite-spin fermions 
residing on nearest-neighbor sites, which corresponds to 
the form (||) of the potential. It is for these reasons that 
we have identified the region in the upper-right corner 
of Figs. 4a and 4b (delimited by a long-dashed line) as 
"correlated bosons" rather than BCS-like, even though 
kF^pair becomes definitely larger than 2ir in this region. 

n 

Another relevant difference of the d- from the s-wave 
solution, is the fact that in the BE region of Figs. 4a and 
4b (— 2/i) does not reduce to the bound-state energy e c of 
the associated two-body problem, unless n = strictly. 
This additional feature of the d-wave solution can also 
be associated with the augmented correlation among the 
composite bosons with increasing n, as discussed previ- 
ously. Such a correlation makes, in fact, the energy re- 
quired to extract two fermions from the system different 
from the energy required to break up a single compos- 
ite boson in isolation, owing to the additional correlation 
energy among the composite bosons. 

The relevance of the correlation energy among the com- 
posite bosons is also suggested by the behavior of the 
"bosonic" chemical potential [Ib — 2/i + e G in the BE 
region (delimited by n < 0.016, irrespective of the single- 
particle dispersion). In this low-density region one finds 
fx B = (U + 2V)n. When < U < 2\V\, \i B is nega- 
tive, corresponding to an effective average attraction be- 
tween the composite bosons. The compressibility of the 
system is also negative at this order, thus indicating a 
tendency toward phase separation. M When U > 2\V\, 
on the other hand, /is and the compressibility are both 
positive; in this case, the effective bosonic average inter- 
action would be repulsive, with an increasing repulsion 
between the composite bosons as \V\ (and U) increases. 
This situation has to be contrasted with the negative-?/ 
Hubbard model, for which the composite bosons become 
asymptotically free as \U\ — > oo. E^j 



IV. CONCLUDING REMARKS 



In this paper we have examined how the particle den- 
sity influences the BCS-BE crossover, by analyzing sev- 
eral types of fermionic interaction potentials both in 
three dimensions (continuum case) and in two dimen- 
sions (continuum and lattice cases). We have reached 
the conclusion that the finite range of the potential as 
well as the absence of a threshold (for a bound state to 
occur in the associated two-body problem) favor quite 
generally the density-induced crossover. 

In particular, for the continuum case we have verified 
that in three dimensions (where a finite threshold exists 
in the two-body problem) it is the finite range of the po- 
tential to make the density-induced crossover possible, 
while in two dimensions the absence of a threshold suf- 
fices to the purpose even when a zero-range potential is 
considered. In addition, for the continuum case we have 
argued that the density-induced crossover (for a finite- 
range potential) occurs at arbitrary values of the inter- 
action strength, owing to the fact that, in this case, the 
particle density can be augmented without bound. 

For the lattice case, we have considered two dimen- 
sions only and investigated, alternatively, the absence or 
presence of a threshold, depending on the symmetry of 
the gap parameter and on the shape of the single-particle 
dispersion relation. We have found that the absence of a 
threshold definitely favors the density-induced BCS-BE 
crossover, and that spreading the range of the poten- 
tial on the lattice makes the density-induced BCS-BE 
crossover to occur over a wider density range. In the 
presence of a threshold (at least for the specific potential 
extending over nearest-neighbor sites that we have con- 
sidered), on the other hand, the density- induced BCS-BE 
crossover turned out not to be possible, owing essentially 
to the absence of the reentrant shape of the curve de- 
limiting the BCS-like region. This situation contrasts 
our finding for the continuum case in the presence of a 
threshold. 

In this respect, an important difference between the 
continuum and lattice cases appears to be the occurrence 
of an intrinsic upper value for the density in the lattice 
case, at least when one considers a simple band only. In 
this case, the Fermi wave vector kp cannot exceed an up- 
per bound of the order it (in units of the inverse of the 
lattice constant). For the nearest-neighbor interaction 
which we have considered in the lattice case, the char- 
acteristic wave vector k , too, can be taken of the order 
7r and the condition kp » k Q (which in the continuum 
case was associated with the reentrant shape) cannot be 
satisfied. Consequently, no reentrant shape of the curve 
delimiting the BCS-like region is expected for the lattice 
case. The occurrence of the reentrant shape, however, 



8 



may not a priori be excluded when considering an in- 
teraction potential extending over distant neighbors in 
the lattice, thus decreasing k a accordingly. In any case, 
the occurrence of an intrinsic upper bound on the den- 
sity (and thus on fcp) will make the BCS-like region to 
disappear for large enough interaction strength, thus pre- 
venting the density-induced BCS-BE crossover. 

We have also found that the presence of an upper 
bound on the density in the lattice case and the finite 
size of the composite bosons for the d-wave solution in- 
duce definite correlations among the composite bosons, 
giving rise to a magnetic correlated state near half fill- 
ing. The simultaneous presence of superconducting off- 
diagonal long-range order and of antiferromagnetic cor- 
relations without long-range order, which we have found 
for the d-wave solution, represents per se an appeal- 
ing result, in the light of the current phenomenology 
of the cuprate superconductors. The interesting physics 
that has emerged for the d-wave solution in the lattice 
case may then provide a general framework for future 
investigations on the BCS-BE crossover, for instance, 
considering finite-temperature effects and going beyond 
the mean-field approximation. In this respect, facing a 
proper treatment of the residual boson-boson interaction 
in the BE region seems to be unavoidable, not only in 
the three-dimensional continuum case (as shown already 
in Ref. fL5fl), but especially in the two-dimensional lat- 
tice case for the d-wave solution, where the boson-boson 
interaction gets considerably enhanced. 

Finally, we point out that the need of a reduced spa- 
tial dimensionality and of a finite range of the potential, 
which we have found for the density-induced BCS-BE 
crossover to occur, matches the generic features observed 
for the cuprate superconductors and strengthen accord- 
ingly the BCS-BE crossover as a possible scenario for 
describing the evolution from overdoped to underdoped 
cuprates. 
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